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Abstract 



C\ \ We compute, from first principles, the quantum fluctuations about instan- 

^ ' ton saddle points of the Euclidean path integral for Einstein gravity coupled 

to a scalar field. The Euclidean two-point correlator is analytically contin- 
ued into the Lorentzian region where it describes the quantum mechanical 
vacuum fluctuations in the state described by no boundary proposal initial 
K^ , conditions. We concentrate on the density perturbations in open inflationary 

ly^ ' universes produced from cosmological instantons, describing the differences 

^^ ■ between non-singular Coleman-De Luccia and singular Hawking-Turok in- 

^\j ■ stantons. We show how the Euclidean path integral uniquely specifies the 

^^ , fluctuations in both cases. 

o\ 

^ ■ I. INTRODUCTION 

o 

H 



The idea that the early universe underwent a period of accelerating expansion i.e. infia- 

^ \ tion, is an attractive one. Such a period of inflation would explain the observed smoothness 

' and flatness of today's universe. But it might also explain the inhomogeneities present in the 

universe. During inflation, quantum mechanical vacuum fluctuations in various flelds would 

^ ■ have been amplifled and stretched to macroscopic length scales, to later seed the growth of 

j^ ■ large scale structures like those we see today. 

However, the theoretical foundations of inflation are still unclear. There is no compelling 
connection between the fleld needed to drive inflation and fundamental theory. The initial 
conditions are usually imposed by hand, or via a somewhat handwaving appeal to primordial 
chaos. The problem of the initial singularity is not addressed. And the usual calculation 
of the quantum mechanical perturbations relies on a rough argument that the perturbation 
modes should be close to the Minkowski space-time ground state when their wavelength is 
far beneath the Hubble radius during inflation [|I|]. This argument is reasonable but hardly 
rigorous since the wavelengths of interest today are typically sub-Planckian during the early 
stages of inflation. 
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An alternative to the usual approach is to confront the problem of initial conditions 
head on by making an ansatz for the initial quantum state of the universe. The Euclidean 
no boundary proposal due to Hartle and Hawking [0 represents one such attempt, and we 
shall implement this proposal in the current work, albeit with slightly different emphasis. 
The main idea is that the path integral may be used to define its own initial conditions, if 
we 'round off' all Lorentzian four-geometries on Euclidean compact four-geometries. This 
is appealing because it is a natural generalization of the imaginary time formalism well 
established in statistical physics, and in some sense represents an unbiased sum over possible 
initial quantum states. One might also add that in field theories and in the theory of random 
geometries in general the only nonperturbative (i.e. lattice) formulations are framed in 
Euclidean terms. The Euclidean approach to quantum gravity is therefore probably the 
most conservative approach to quantum gravity, building upon techniques which are well 
proven in other fields. 

Recently, Hawking and one of us found a class of singular but finite-action Euclidean 
instantons which can be used to define the initial conditions for realistic inflationary universes 
0] . For a generic inflationary scalar potential there exists a one-parameter family of singular 
but finite-action instantons, each allowing an analytic continuation to a real open Lorentzian 
universe. If these instantons are allowed, then open inflation occurs generically and not just 
for potentials with contrived false vacua as was previously believed |^. In this paper, we 
investigate the spectrum of perturbations about such singular instantons as well as the more 
conventional non-singular Coleman-De Luccia instantons |^ previously used to describe 
open inflation. We shall be particularly interested in determining whether there are any 
observable differences between the perturbation spectra produced in these two cases. 

The instanton solution provides the classical background with respect to which the quan- 
tum fluctuations are defined. In the Euclidean path integral approach, one can in principle 
compute correlators of the quantum fluctuations perturbatively to any desired order in h. 
Unlike the usual approach to inflation, there is no ambiguity in the choice of initial conditions 
at all, because the Euclidean Green function is unique. We shall compute the Euclidean 
Green function here for instantons of the singular type, as well as for regular Coleman-De 
Luccia instantons such as occur in theories with a false vacuum. 

The question of whether singular instantons are allowed has provoked some debate in 
the literature (see 0, 101! and references therein). In this paper we add to the evidence 
in their favour by showing that the spectrum of fluctuations is well defined in the presence 
of the singularity. The Euclidean action by itself specifies the allowed perturbation modes 
without any extra input. In a parallel paper it is shown that the same is true for tensor 
perturbations. These calculations demonstrate that at least to first order in h the quantum 
fluctuations about singular instantons appear healthy. In principle, the present framework 
also offers a method for computing higher order corrections, although one expects that 
the non-renormalizability of gravity would introduce new free parameters describing the 
coefficients of the higher order counterterms. 

There already exists an extensive literature on the problem of fluctuations in open in- 
flation. There have been a number of difficulties in obtaining precise results and to our 
knowledge there are no calculations yet available for generic potentials which do not make 
one approximation or another (see e.g. [|1^, [|TT|, |]12|)- There are also a number of unre- 



solved ambiguities which arise in precisely which modes to allow (see the discussion in |]T^ 



and references therein). 

We believe we have isolated and overcome these problems in the present work and re- 
duced the problem to a complete prescription which may be numerically implemented. The 
main idea of our method is to compute real-space correlators in the Euclidean region and 
analytically continue them to the Lorentzian region. In contrast, previous work has treated 
the perturbation modes of each wavenumber separately. This can lead to confusion since 
the naive open universe modes provide only a partial description (the 'sub-curvature' piece) 
of the relevant correlators. The remaining 'super- curvature' piece is not expressible in terms 
of these modes. In our method, both pieces are automatically included, and the connection 
between them is thereby clarified. 

In the Euclidean real-space approach the problem of sub-Planckian modes does not occur 
as directly as in the usual approach, since the initial conditions are defined in a manner that 
makes no mention of the mode decomposition. In analogy with black hole physics, there 
is reason to hope that the results obtained are likely to be insensitive to short distance 
sub-Planckian physics. 

Euclidean quantum gravity is not a complete theory. As well as being non-renormalizable, 
it suffers from the well known conformal factor problem, the Euclidean Einstein action being 
unbounded below. This problem does not affect the calculations reported here because the 
spatially inhomogeneous perturbations have positive Euclidean action. It is only the spatially 
homogeneous modes which suffer from the conformal factor problem. Even for these modes 
we think that there are no grounds for pessimism. For pure gravity, or for gravity with 
a cosmological constant, the conformal factor problem disappears when the gauge fixing 
procedures are carefully followed. The technology for scalar fields coupled to gravity has not 
yet been worked out, but as far as we know there is no insuperable obstacle to doing so. 

Although the simplest (e.g. 0^ or 0'^) inflationary potentials yield in this approach 0] a 
most probable universe which is much too open to be compatible with observation, there are 
other scalar potentials which yield acceptable values closer to unity [0]. In this paper we 
do not address this question, but formulate our results so that they apply for an arbitrary 
scalar potential. 

Finally, we note that in a generic inflationary theory there also exist singular instan- 
tons yielding real Lorentzian closed inflating universes. We shall investigate the quantum 
fluctuations about such instantons in a future publication. 

II. FLUCTUATIONS FROM THE EUCLIDEAN PATH INTEGRAL 

In this section we discuss the principles of our method, postponing the technicalities to 
the next section. We choose to frame the Euclidean no boundary proposal in the following 
form. We take it to be an essentially topological prescription for the lower limit of the 
functional integral. We write the quantum mechanical amplitude for the state described by 
three-metric g^ and matter fields </> as 

W.^]-] [Dg]W]e^^ (1) 

and our prescription is to integrate over all Euclidean/Lorentzian four-geometries g of the 
form shown in Figure ^ with associated matter fields 0, bounded by the three-geometry 




FIG. 1. The quantum amplitude for a three-geometry Sj is given as the path integral over all 
Lorentzian geometries matched on a three-surface S to a compact Euclidean four-manifold. All 
geometries of this type are to be integrated over. 

Sj and matter fields present in the final state. The Euclidean region is essential because 
there is no way to 'round off' a Lorentzian manifold without introducing a boundary. We 
sum over all matching surfaces S and Euclidean regions bounded by S. Our prescription 
differs from that of Hartle and Hawking f^ in that we shall not impose regularity of the 
four-geometries summed over. Such a prescription seems to us to be at odds with the basic 
principles of path integration. Regularity of the background geometry and the fluctuation 
modes should emerge as an output of the path integral rather than be input, and we shall 
see this occur in our calculations below. 

The expression (0) is only formal, and part of our investigation will be to see whether we 
can calculate it in a perturbative expansion around saddle point solutions i.e. instantons. 
The relevant instantons are real solutions of the Euclidean field equations which possess a 
surface E on which they may be analytically continued to a real Lorentzian spacetime. The 
condition for this to be possible is that the normal derivatives of the three-metric and matter 
fields should be zero on S. Note that regularity of the four-geometry of the instanton, needed 
here so that analytic continuation is possible, is a consequence of the instanton being a saddle 
point in the sum over all four-geometries and therefore a solution of a partial differential 
equation. Such analytic four-geometries are of course a set of measure zero in the original 
path integral. 

We wish to calculate correlators of physical observables in the Lorentzian region. Each 
instanton provides a zeroth order approximation, giving us a classical background within 
which quantum fluctuations propagate. To first order in h the quantum fluctuations are 
specified by a Gaussian integral. One can perform the integral in the Euclidean region 
and then analytically continue in the coordinates of the background solution to find the 
quantum correlators in the Lorentzian region. As we shall now see, the analytic continuation 
is guaranteed to give real-valued Lorentzian fields and momenta if the background solution 
is real in both the Euclidean and Lorentzian regions. 

The quantum mechanical amplitude for the fluctuations about a particular background 
solution B is given from (|l|) by (henceforth we set ^ = 1) 



^s[^^ 0] ~ e'^^^^^'*^^ f ^ ' [VSg] [V6^] e'^^^^>^'^*^ (2) 

where the metric g = gs + ^g and fields 4> = (pB + ^(p- <S2 is the action for second and higher 
order fiuctuations. Let P(xi) and Q(x2) be two observables at xi and X2 on E/. Then their 
correlator is given by integrating |\l/p times V{xi)Q{x2) over a complete set of observables 
0{x) on a final state three-surface Ej containing xi and X2, 

(P(a;i)Q(a;2)) = U j dB j [VO{T.f)] ^>b [O]* ^> b [O] V{x^)Q{x2) (3) 

with M an appropriate normalization constant. The full quantum correlator involves the 
sum over all background solutions B. If no other constraints are imposed, the solution with 
the lowest Euclidean action will dominate. 

The insertion of the sum over a complete set of states means that the result @ may be 
viewed a single 'doubled' path integral. The 'initial' state is established on the Euclidean 
region bounded by the first version of S. There is then a Lorentzian region running forward 
with weighting factor e**^ to Ej on which the operators of interest are located. Then there 
is a region running back to a second version of S, with weighting factor e"**^, and finally 
the geometry closes on another Euclidean compact four-manifold. In the instanton approx- 
imation, the two Euclidean regions are actually copies of the same half-instanton. But the 
quantum fiuctuations on the upper and lower halves are independent - the wavefunction \& 
and its conjugate \1'* involve independent integrations over the perturbations. 

In formula @ we can now continue Xi and X2 back into the Euclidean region of the 
background solution. If the Lorentzian continuation of the instanton is real, then as men- 
tioned the Lorentzian part of the path integral involves a factor of e*'^ coming from \1/ and 
a factor of e~^^ from \1/*. These cancel exactly allowing us to deform Ej back towards the 
Euclidean region. Upon reaching the Euclidean region however, we discover that both \l/ 
and \1/* involve e""^^ so there is no cancellation. In the correlator we are left with an overall 
e"^"^^ where Se is the half-instanton action, and the action for the fiuctuations is just the 
Euclidean action evaluated over the doubled half-instanton. Correlators calculated in the 
Euclidean region and continued to the Lorentzian region are therefore equal to those com- 
puted from the Lorentzian path integral with Euclidean no boundary initial conditions only 
if the background solution is real. Note that the cancellation of e*"^ and e"*"^ occurs for any 
three-surface Sj in the Lorentzian region. There is no requirement that Sj be spacelike. In 
fact, there is no reason for a complete surface Sj to exist at all. If one computes correlators 
in the Euclidean region as we shall, all that matters is that a smooth continuation of local 
observables exists into the Lorentzian region. If there are singularities, they can be avoided 
by choosing an appropriate continuation contour. 

The above construction is closely analogous to the imaginary time formalism for ther- 
mal field theory, where real-time correlators can be calculated by analytic continuation of 
Euclidean ones. In the present context, the compact nature of the Euclidean instantons 
has the same effect on correlators that the periodicity in imaginary time does in thermal 
field theory. The periodicity of the instanton solution introduces thermal weighting factors 
corresponding to the Hawking temperature of the background spacetime. 




FIG. 2. The Euclidean/Lorentzian manifold in the vicinity of the regular pole. 
III. INSTANTONS AND ANALYTIC CONTINUATION 



Let us briefly recall the form of the instantons we are interested in. We consider Einstein 
gravity coupled to a single scalar field with potential energy V{(j)). We seek finite-action 
solutions of the Euclidean equations of motion. If V{(j)) has a positive extremum, there is 
a solution in the form of a four-sphere. This has the maximal symmetry allowed in four 
dimensions, namely 0(5). In general however, no such instanton exists, and the highest 
symmetry possible is 0(4). The instantons are described by the line element 



rfs^ = da^ + b'^{a)dnl = da"^ + 6^ (a) (rf^^ + sm^{'4))dnl 



(4) 



with h{a) the radius of the three-sphere. The Einstein and scalar field equations take the 
form 



-^' 



Vj 



K,, 



V 



(5) 



where n = SnG, (p^a = da(p and V^j, = ^. 

The point of |^ was that for a gently sloping potential of the kind needed for inflation, 
there exists a one-parameter family of solutions to these equations labelled by the value of the 
scalar field at the north pole of the instanton. At the north pole the metric and scalar field 
are regular. At the south pole, which is singular, the scale factor goes to zero as (a^ — o") ^ , 
and the scalar field diverges logarithmically in a. These singular instantons are in general 
not stationary points of the action and they should be interpreted as constrained instantons, 
where the constraint is imposed on a small three-surface surrounding the singularity |§]. 

The various analytic continuations are fixed by the fact that the Euclidean/Lorentzian 
geometry in the neighbourhood of the north pole is locally flat. The coordinates we use for 
the instanton and its continuation reduce in this neighbourhood to those used in mapping 
flat space onto the Milne universe, see Figure ^ There is then an obvious choice of regular 
coordinates which allows one to uniquely fix the required continuations. 

Consider flat Minkowski space with line element —dT"^ + dR'^ + R^dVL\. The analytic 
continuation to Euclidean time is performed by setting Te = iT . Under this continuation 



the Lorentzian action S yields the positive Euchdean action, iS = —Se- We may describe the 
situation in 0(4) invariant coordinates Te = acosr^, R = asinQ, with Q the polar angle 
on the three-sphere. These coordinates correspond to those used for our 0(4) invariant 
instantons. As Q runs from vr to zero, Te runs from —a to a. For fixed a, the Euclidean 
action involves the integral /^g. (ITe- The continuation to the Lorentzian region is described 
by distorting the Te contour into the complex plane. The new contour runs from —a to zero, 
up the positive imaginary axis where the e""^-^ factor gives e*"^ needed for the wavefunction \E' 
in (^. The continuation then runs back down the imaginary axis, giving e~*'^ and along the 
real axis from zero to +cr. If we translate Te into Q, the first two segments of the contour 
run from Q = tt to Q = 7i/2 and then down the line Q = {tt/2) — it', with t' real and positive. 
The latter segment gives us the e*"^ term in the 'doubled' path integral. 

Following this continuation on the surface Te = 0, from the Euclidean region we obtain 
a Lorentzian spacetime with timelike coordinate t'. We have T = cr sinht' and R = a cosh t'. 
Since these formulae imply R^ — T^ > 0, these coordinates cover only part of the Lorentzian 
region, the exterior of the light cone emanating from the origin of spherical coordinates at 
T = 0. We call this region II. We can then perform a second continuation t' = x — z| by 
setting T = tcoshx and R = tsinhx to obtain the metric in the interior of the light cone, 
which we call region I. In the vicinity of the regular pole, these new coordinates are those 
of the Milne universe. Globally, region I is the inflating open universe. 

Hence the required continuations are 

Euclidean ^ I : a = it, Vl = —ix, 

Euclidean ^ II : Q = -- it' , (6) 

yielding the following line elements 

Euclidean : dcr'^ + b'^{a)dQl, 

Region I : - dt^ + a? (t) {dx^ + sinh^ xd^l), 

Region II : da^ + b\a)l-dt'^ + cosh^ t'dnj). (7) 

The function 6(cr) ^ a as a ^ 0, and a{t) -^ t as t — > 0, so it follows from (|]) that b = ia. 

The Euclidean action is expressed as an integral over the coordinates a and Vt, and 
then regarded as a four dimensional contour integral. We may then distort the integration 
contour into the regions of the complex coordinate space corresponding to the Lorentzian 
regions I and II. This procedure uniquely defines the path integral measure for Lorentzian 
correlators. 

In what follows, it will be very useful to work in terms of a conformal spatial coordinate 
in the Euclidean region and in region II, defined by 

For singular instantons, X = corresponds to the singular pole and X — > cxd corresponds 
to the regular pole. For regular (Coleman-De Luccia) instantons, the second regular pole is 
at X ^ — cxD. 

It is also useful to work in a conformal time coordinate r in the open universe, defined 
to obey dt = a{t)dT. For many purposes we shall need to relate r to X. To do this we 
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FIG. 3. Contour for X (t). 

extend our integral definition of X in equation (|^) into tlie complex a-plane. The required 
integration contour is shown in Figure ^. One has 



where we define the lower limit of the conformal time r via 

- da /•* dt' \ 




b{a) Je a{t' 



(10) 



so that r runs from — oo at the beginning of the open universe to an approximate constant 
when inflation is well underway, and then finally diverges logarithmically to cxd at late times 
in the open universe when a (t) ~ t. 

The conformal structure of the Lorentzian region is shown in Figure ^. As the diagram 
indicates, the singularity is timelike and visible from within the open universe region, labelled 
I. It is interesting to ask when the singularity first becomes visible to an observer in region 
I. Any such observer can by symmetry be placed at the origin of spherical coordinates. To 
find the null geodesies, we first change coordinates to conformal time in region I, defined in 
equation (|T^. Null geodesies incident on the origin of spherical coordinates, x = 0, at a 
conformal time Tq obey 

X = ro - r. (11) 

With the above continuations @ the conformal space coordinate X in region II obeys 

X = t'- To. (12) 

The singularity is located at X = 0, and the transition from the Lorentzian to the Euclidean 
region is at t' = 0. It follows that singularity first becomes visible in the open universe when 
tq = 0. As mentioned above, with the above definition of r, infiation ends at negative 
conformal time. For example, for a ^ni?(fp' potential with 60 efolds of infiation, we find 
the end of infiation occurs at r = —1.70 in units where the space curvature is minus one. 
Evolving forward into the late universe, one finds that the singularity is first visible at late 
times when universe is becoming curvature dominated. 

Although the singularity is visible within the open universe, and has a definite observable 
effect, as we shall see the quantum fluctuations are nevertheless well deflned in its presence. 
The singularity acts as a reflecting boundary [|15[, and nothing enters the universe from it. 

8 




FIG. 4. Conformal spacetime structure of the classical background solution. 

IV. THE PATH INTEGRAL FOR SCALAR PERTURBATIONS 

In this section we derive the action appropriate for fluctuations about the instantons 
described above. Since the background solution satisfles the fleld equations (including an 
appropriate constraint if that is required ^), the leading term in the action occurs at second 
order in the perturbations. We shall calculate this second order term and the perform the 
Gaussian path integral determining the quantum correlators to flrst order in h. 

We compute the relevant action in the open universe region, where a Hamiltonian treat- 
ment is straightforward. Our discussion follows the notation of ref. [|I^, although we shall 
work strictly from the path integral point of view. The action for scalar perturbations re- 
duces to that for a single gauge-invariant variable q, related to the Newtonian potential "^n 
standardly used in the analysis of inflationary quantum fluctuations. We analytically con- 
tinue the variable q into the Euclidean region of the instanton, where its action is positive. 
The real-space correlator of q is computed in the Euclidean region and flnally analytically 
continued back into the Lorentzian region. 

Following standard notation we write the perturbed line element and scalar fleld as 



ds^ = a^ (r) ( - (1 + 2 A) dr^ + Sidx'dr + {%j + hij) dx'dx^ 



(13) 



where jij is the metric on the background three-space, r is the conformal time, 0o (t) is the 
background acalar fleld and a{T) is the background scale factor. 



We decompose Si and hij as follows (see |T^) 



h 



*j 



-2ip-fij + 2E\ij + 2F(i\j) + t 



t]i 



Si = B\i + Vi. 



(14) 



Here | denotes covariant derivative on the background three-space, if), B and E are 
scalars, Vi and Fi are divergenceless vectors, and tij is a transverse traceless tensor. In 
general, with a suitable asymptotic decay condition, the above decomposition is unique up 
to B ^ fi+constant, Fi -^ Fi + Ki, with Ki a Killing vector. For compact 3-spaces there 



9 



is in addition an ambiguity since the metric is unchanged under ip^ip — (,E^E + (, 
with (A + 3) C = 0. However for the compact Euchdean instantons we consider all modes 
obeying (A + 3) C = are actually pure gauge and we shall in any case have to project them 
out. 

Under an infinitesimal scalar coordinate transformation x^ — > x^ + X'^, where A'^ = 
(A", A'M, the perturbations in (p!4D transform as 

ip^ij- 7^A°, 5 ^ S + A' - A°, 
A-^A + A°' + 7^A°, E^E + X, (50 -^ (50 + 0[)A° (15) 

where Ti = a' /a, and here and below prime denotes derivative with respect to conformal time. 
The action is invariant under these transformations. We must pick a gauge in order to fix 
this invariance and obtain a unique result. As is well known, the computation of infiationary 
perturbations from a single scalar field is simplest in conformal Newtonian gauge, defined by 
setting B = E = 0. This completely fixes the gauge freedom. Equivalently, one can define 
the gauge-invariant variable [^ 

xij^ = ^lj-n{B-E'). (16) 

As long as there are no anisotropic stresses, "^^ is governed by a second order differential 
equation in time, and all perturbations are determined from \E'7v with the use of the Einstein 
constraint equations. An approximately-conserved quantity x can be constructed from from 
\I/7V which may be used to match the super- Hubble radius perturbations across the reheating 
surface and into the late universe relevant for observation (see e.g. [^). 

In this section we derive the path integral appropriate for computing "^^ correlators in 
the open universe. We shall do so in a manifestly gauge-invariant manner. The first point 
to note is that "^n involves a field velocity and in the path integral formalism the first step 
is to convert this to a canonical momentum. This requires us to use the Hamiltonian (first 
order) form of the path integral. A second merit of this formalism is that integration over 
the nondynamical lapse and shift fields imposes the Einstein constraint equations as delta 
functionals, enabling further integrations to be performed. Our discussion parallels that 



of Appendix B of reference |jT3|, but is slightly more concise. We shall also be careful to 
keep certain surface terms which determine the allowed fluctuation modes about singular 
instantons. 

Our starting point is the action for gravity plus a scalar field 

5 = ^ / d'xy^ (i? - i V^0V'^0 - V (0)) - ^ I £x^K, (17) 

where K is the trace of the extrinsic curvature of the boundary three-surface. The surface 
term is needed to remove second derivatives from the action, so that fiuctuation variables 
are constrained on the boundary but their derivatives are not. The decomposition ([M]) 
is substituted into equation ([17|) , keeping all terms to second order. The scalar, vector 
and tensor components decouple. The vector perturbations are uninteresting because the 
Einstein constraints force them to be zero. The tensor perturbations are discussed in the 
parallel paper M. The second order action for the scalar perturbations reads 
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S2 = yJ drd^xa^^l - 6^'^ - UHA^' + 2AV^ {2 A - ^) - 2 (W + 271^) A^ 

+/C (-6^2 ^ 2^2 + 12^A + 2 (5 - E') A (5 - F)) 

+4A (5 - E') (|0'o^0 - V^' - HA)} (18) 

where A is the Laplacian on the three-space. Integrations by parts in the spatial directions 
may be freely used in anticipation of the fact that the fluctuations are determined in the 
Euclidean region where the three-space is an 5*^ so there are no surface terms. We shall 
eventually need to be more careful about integrating by parts with respect to r, at least in 
the case of singular instantons, because r continues to the coordinate X in the Euclidean 
region which terminates on the singular boundary. However, even here we may integrate by 
parts freely until we have expressed our action density in Hamiltonian form for the observable 
of interest, namely \Ef at. After this last stage we need to be careful to retain all surface terms. 
As we shall see, after continuation to the Euclidean region these surface terms determine 
the set of allowed fluctuation modes. 

We introduce the momenta canonically conjugate to ip, E, and Scj) as 

n^ = '^iVl f^^' + 3-^'^S(p - SUA - A (5 - E')\ , 

Ue = ^^-^ (ij' - ^^'o5^ + nA- JCiB - E')y 

Ils^ = a^Vli^<f>'-<P'oA)- (19) 

We shall only consider modes for which —A — 3/C is positive, and in this case 11^ and 11^ 
are independent and we can solve for the velocities in terms of fields and momenta. As an 
aside, we mention a subtlety associated with the modes known in the literature as 'bubble 
wall' fluctuation modes [jl^, [0. These modes have have A + 3/C = 0. Such perturbation 



modes are possible in the open universe in spite of the fact that they possess the 'wrong sign' 
for the Laplacian and therefore grow exponentially with comoving radius. However, if one 
expands in harmonics, for / < 2 such modes may be gauged away. And modes with / > 2 
are singular when continued into the Euclidean region, since with A^; = — (n^ — 1) = —3 
the regular eigenf unctions all have I < n = 2. So there are no physical fluctuation modes 
with A + 3/C = 0. However, our gauge-invariant action will be zero for the A -|- 3/C = 
gauge modes, and we will need to project out their contribution at various stages of the 
calculation. 
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Proceeding with this caveat, we rewrite the action ([T8| ) in first order form 



S2 = J dTd'xlu^^' + UeE' + Us^6(f)' 



K ( 2 oTT TT 3^2 2(A + 3/C)^2 



4a2 V7 (A + 3^) V"^"^ ^ ^"'"'' ^ A^^ ^ 



K 



n^^ J - 2 '^on^.^'^ 



aV7 



^ ("(A + 3^) ij^-'^U + ?,ic-n^-n' + ^\ 6 A 



-BUe - A (-mi^ + 0'on5<^ + ^^^ (^_ ( A + 3/C) V^ + I (7^0[, - 0'o') '^<^) ) I • (20) 

As stated above, we would hke to evaluate the correlator for the gauge-invariant variable 
\E'7v. We may express \E'7v as a dynamical variable in terms of fields and momenta as follows 

Tin /An^ + 3n£\ ,^. 

This is singular for A = so our discussion will only be valid for the inhomogeneous 
perturbations, which is all we shall calculate here. The homogeneous perturbation modes 
require a separate treatment because the Euclidean action is not positive and to cure this 
the conformal factor must be decoupled from the physical degrees of freedom. We do not 
address this problem here, but shall simply project out the homogeneous modes on the S"^ 
or H^ at each stage of the calculation. 

We now add a source term —i J Jn'^n to the quadratic action, and perform the path 
integral over all fields and momenta. The non-dynamical variables A and B only occur 
linearly in 5*2, and functional integration over these fields gives us delta functionals imposing 
the Goi and Gqo Einstein constraints. We use these delta functionals to perform the He and 
n^,^ integrals. This sets He to zero in the term multiplying the source. There is then no 
residual E dependence in the action and and the functional integral over E is ignored as an 
infinite gauge orbit volume. 

There is still a residual gauge freedom in the action, corresponding to A° reparametriza- 
tions (see above). However the action density can be expressed, up to a total derivative, in 
terms of "^n and its canonically conjugate gauge-invariant variable 

Hjv = -n^ - -^ (A + 3/C) 7/. ^ (A + 3/C) (50. (22) 

The action density is now independent of 6(j), so this too can be integrated out as a gauge 
orbit volume. 
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It is convenient to rescale the coordinate a = — S-\l/ 
the source Jn = ^Jh 



AT, the momentum p = -^IItv, and 



,q — 2^ ^N- We then perform an integration by parts to represent the action in 
canonical form S = J [pq' — H{p, q)]. From now on, all surface terms must be kept. Finally, 
we perform the Gaussian integral over the momentum variable p to obtain the reduced action 
for q, 

iS2 (J) =-\j dTd?x^l (A + 3/C) g (6 + A + 3/C) q 



+ (A + 3X:) 



qq + 




+ 



jj,q 



where 



O 




(23) 



(24) 



The result for the bulk term agrees, up to a sign, with that obtained in |T^. We are now 
ready to analytically continue to the Euclidean region and perform the Gaussian integral to 
obtain the qq correlator. 



V. THE EUCLIDEAN GREEN FUNCTION 

The continuation to the Euclidean region is performed as described above, in equation 
(^). The Laplacian A continues to —A3, where A3 is the Laplacian on S"^ and the constant 
/C continues to itself. We set /C = — 1 from now on. Finally, the variable q continues to 
itself. 

The Euclidean action is given by 



^^2 (J) = -Sl^ (J) = -- / dXd'n^{ (-A3 - 3) g (O - A3 - 3) g 



+ (-A3 - 3) 






JqQ 



(25) 



where the volume element on the three-sphere is d^Q^ and now 

CI rL,,/>o ,./X\ Qj 



O 




rfX2 



+ U{X). 



(26) 



where primes now denote derivatives with respect to X. The bulk term in the action is 
positive definite for the inhomogeneous modes of interest, as was noted by Lavrelashvili 

The path integral over q is Gaussian and is performed by solving the classical field 
equation for q{J) and substituting back. After the appropriate normalization, we obtain the 
generating functional for q correlators, exp(| // JGeJ), where Ge is the Euclidean Green 
function. The two point correlator is then {q (X) q (X')) = Ge{^, X'). The Euclidean Green 
function is the solution to the equation 
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(-As - 3) (6 - A3 - 3) Ge = g-h^'\X - X') - <5(X - X')^^. (27) 

Here X = (X, x*), and g is the determinant of the metric. Q is the angle between the two 
points X and x' on the three-sphere. The second term on the right hand side of equation (pTl) 
is present because we project out the —A3 = 3 gauge modes as mentioned above. Its 
coefficient is fixed by the requirement that the right hand side be orthogonal to the —A3 = 
3 modes under integration over the three-sphere. We shall also project out the A = 
homogeneous mode, but for the sake of brevity we shall not write out the relevant terms 
explicitly. 

Equation ( pT]) is solved by expressing both 6{X — X') and Ge as sums over a complete 
set of eigenmodes of O and equating coefficients. O is a Schrodinger operator, and its 
eigenfunctions obey 

OMX) = (-^ + f/(x)) MX) = (/ + ^)Mx), (28) 

where the potential U{X) is given in equation (^6]). Near the regular pole of the instanton 
i.e. X -^ 00, we have 0o,o- ~ cr ~ g~^ , so that 0q ~ e~^"^, and U ^ 4. There is therefore 
a positive continuum starting at p^ = 0. For gently sloping inflationary potentials, there is 
also generally a single bound state with — 1 < p^ < 0. 

For singular instantons, near X = the potential term is repulsive and diverges as ^^, 
as noted by Garriga [|l^. There are two sets of eigenmodes of O for each p, behaving as 
X~i/2 or as X^/^ for small X respectively. Now the importance of keeping the surface terms 
in the Euclidean action (^Sf) becomes clear. The surface terms are positive inflnity for the 
divergent modes. Thus the Euclidean action by itself completely determines the allowed 
spectrum of fluctuations. Regularity of the eigenmodes does not need to be imposed as an 
additional, external condition. Only fluctuations vanishing at the singularity are allowed, so 
in effect the singularity enforces Dirichlet boundary conditions. 

Now consider for fixed real p the solution to the Schrodinger equation (|28| ) which behaves 
as X2 for X — >■ 0, which we shall denote il)p{X). Being the solution to a differential equation 
with finite coefficients, this is analytic for all finite p in the complex p-plane |^ . It is also 



useful to define the Jost function gp{X) as the solution which tends to e"^'^^ as X ^ 00. This 



is analytic in the upper half complex p-plane, as seen by iterating the integral equation |21 



1 /""^ 
gp{X) = e^P^ + - / sinpiY - X) V{Y) gp{Y) dV. (29) 



p JX 

The two solutions are related via 



MX) = a^gpiX) + bpg^X). (30) 

Since the differential equation is even in p, ip_p{X) = ippi^X) and bp = a-p. Completeness of 
the eigenfunctions then allows us to write the following representation of the delta function, 



5 (X - X') = r MpM^dp + < (X) < (X') , (31) 
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where we have included the contribution of a single assumed normalized bound state wave- 
function iIj^{X). One may check the normalization of the continuum contribution in equa- 
tion (^) as follows. For X and X' large, substituting in (|30| ) one sees that the terms going 
as e^^pi^-^') integrate to give the correctly normalized delta function. 

It is also instructive to see how the right hand side vanishes for any X ^ X' . For 
definiteness let us take X > X'. Substituting the decomposition (^) into the integral for 
ipp{X), we have 



oo Anapa^p J-oo Airapa^p 

gp{X)^p{X' 



27ia^p 



dp. (32) 



We now distort the integral onto a semicircle at infinity in the upper half p-plane. The 
contour at infinity gives no contribution since gp{X)-ipp{X') decays exponentially, as e^^^^"^^ \ 
Since ^p(X') is analytic for all p, and gp{X) is analytic in the upper half p-plane, the only 
contribution to the integral comes from zeroes of a^p in the upper half p-plane. These zeroes 
correspond to bound states as may be seen from the expression (|30|). If a_p vanishes, the 
corresponding solution decays exponentially at large X. Thus for each zero of a^p we have 
a normalizable bound state with eigenvalue p^. But the Schrodinger equation only has real 
eigenvalues, so this is only possible for imaginary p = iA, with A real. 

Next we compute the residue of the pole, which we assume to be simple. From the 
Schrodinger equation one obtains the following identity 

We integrate both sides over X from zero to infinity. The left hand side yields a difference 
of surface terms. Only the surface term at infinity is non-zero, and we evaluate it using 
the asympotic form for tpp. Near p = iA, a^p vanishes as {p — iA) .dpa-p\p=i\. Then the 
integrated equation (B^) gives us 9„a_pL=jA = — — , where the normahzation constant 
M = /o°° iplj^dX. So the contribution from the zero of a_p to the integral (|3^) is just 

2 

-^gi!^{X)gi!^{X') . But ^=gi\{X) is precisely the normalized bound state wavefunction 

ip^i^{X), so the contribution from the zero of a-p is exactly cancelled by the bound state 
term in (|3TD. Hence the right hand side is zero as required. 

All of this is simply understood. Imagine deforming our potential U{X) to one possessing 
no bound state. Then the integral in (|32D, with p running along the real axis, would give 
the delta function with no bound state contribution being needed. Now as one deforms the 
potential back to U{X), a zero in a_p corresponding a bound state would appear first at 
p = 0. This yields a pole in the integrand of (|32|) but the contour may be deformed above 
it since the rest of the integrand is analytic in the upper half p-plane. As the potential 
becomes more and more negative, the zero of a^p moves up the imaginary p axis to p = iA. 
But the integral still equals 5{X — X') as long as one takes the contour above the pole (see 
Figure ^. If one deforms the contour back to the real p axis, the bound state contribution 
discussed above is produced from the residue of the pole. 
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FIG. 5. Contour of integration avoiding the bound state pole. 

In the case of non-singular instantons, a similar procedure may be followed. Here, since 
b{a) goes linearly to zero at each pole, X should be defined as /^* ^t^, where at is say the 
value of a for which b{a) is a maximum. Then X ranges from — cxo to oo, and for each 
value of p^ we need two linearly independent mode functions. These may be taken to be 
gf^ (X), defined to tend to e'^P^ as X ^ -oo, and gf^^^ (X), defined to tend to e'P^ as 
X — > oo. These can be shown to be orthogonal, and analytic in the upper half p-plane. As 
X ^ oo, if we write g^f' (X) -^ CpC'P^ + dpC-'P^ , then as X ^ -oo, g'^'^^^ (X) -^ dpC'P^ - 

c_pe-'P^. Finally, we may express 6 {X - X') as /^ ^^— ^Sf^^P + '^S. (^) ^iX (^') in 
close analogy with the singular case. 

Returning to our discussion of the singular case, from the representation (0) of the delta 
function we may now construct the following ansatz for the Green function 

Ge = r Cp (Q) MpM^dp + Qa {Q) ^^a (X) ^il (X') (34) 



-p^-p 



where as above ^2 is the angle between the two points x* and x* on the three-sphere. 

Substituting (0) into (|27| ) we obtain an equation governing the 'universal' part of the 
correlator, Cp (Q) . This reads 

sin 20 

(A + 3)(A - 1 - p')Cp (O) = 5'in) - ^^, (35) 

where A = 9^ + 2 cot Qdn- 

We first find the particular integral for the term on the far right and then the four linearly 
independent solutions of the homogeneous equation. The latter are 

cos 20 sin 20 coshpO sinhpO , . 

—FT^^Tr^ . , and . (36 

smO smO smO smO 

We combine the 2 solutions singular at O = to cancel their leading singularities, and 
arrange that the term linear in O at small O has the correct coefficient to match the delta 
function in (|35D . Then we choose the coefficient of sinhpO/ sin O to make Cp finite at O = tt 
(i.e. for antipodal points on the 5*^). Finally we choose the coefficient of the sin 20/ sin O 
term to make the projection of Cp onto the gauge modes zero. Hence we find 

1 sinhp(0 — tt) (tt — 0)cos20 



Cpin) 



An [p^ + 4) sinh pvr sin O 47?^ {p^ + 4) sin O 

/I 1 \ sin 20 

V7rMp' + 4)' ~ 167r2(p2 + 4)j sinO ' 
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(37) 



where the first term gives us the contribution one might expect by analogy with the usual 
Euclidean scalar field vacuum. We shall be interested in the behaviour of Cp {VL) in the 
complex p-plane. Here, the role of the extra terms is just to remove the double pole at 
p = 2i. Had we written out the terms involving the homogeneous mode on the three-sphere, 
we would find that their effect would be to remove the pole aX p = i. 

We now have a convergent expression for Ge in the Euclidean region. Note that whilst 
Ge is perhaps most naturally expressed as a sum of regular eigenmodes, discrete because 
the space is compact, we have instead expressed it as an integral. The integral formula is 
more useful for analytic continuation since it is already close to an expression of the form 
we desire in the open universe, namely a integral over a continuum of modes. 

Nevertheless it is interesting to see how the Euclidean Green function appears as a 
discrete sum. For X > X', we may close the p contour above to obtain an infinite sum, 
convergent in the Euclidean region, 

„t^3 47r2 (n2 - 4) sinfi a_„, ' ^ ^ 

where gni{X) ~ e~"^ at large X. This demonstrates that our Green function is analytic in 
proper distance a ~ e~^ at the north pole of the instanton, as it should be. For non-singular 
instantons the argument generalizes to show that the Green function is analytic at the south 
pole too. 



VI. THE LORENTZIAN GREEN FUNCTION 

We now wish to continue our integral formula for the Green function given by (|3^), ( pTf ) 
into the open universe region. This involves setting Q = —ix and continuing the conformal 
coordinate as described in equation (^. 

To perform the continuation we take X > X' and write Ge as 

c 27ra_p 

where the contour C for the p integral has been deformed above the bound state pole as 
described above (see Figure]^). 

We can perform the analytic continuation Q = —ix immediately. However, the X 
continuation is more subtle, because gp ~ e*^^ at large X, and unless we are careful terms 
like ~ e^'^ occur which cause the p integral to diverge. We circumvent this problem as 
follows. For X — X' > we have 

gpiX)^PpiX') e^P^ ^ e~^^ g,,iX)MX') ,^^. 

c 27ra_p (p2 + 4) IQir^ a.2i 

since the integrand is analytic in the upper half p-plane and the integral may be closed 
above. By inserting sinhpvr/sinhpTr = 1 under the integral one sees that the integral (^Of ) 
with a factor e^'^/sinhpvr inserted is equal to that with a factor e~^'^/sinhp7r inserted, plus 
the remaining term on the right hand side. We use this identity to re-express the term from 
(^) behaving as e^'^"'"*^. 
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The term on the right hand side of (^OD may be then combined with the analytic con- 
tinuation of the term involving vr cos 2^2 on the right hand side of (^) to produce a term 
proportional to sinh 2%/ sinh x- The remaining terms in the correlator arising from the sec- 
ond and third terms in (|37|) are proportional to x cosh 2%/ sinh x and sinh 2%/ sinh x- These, 
and the constant terms we have not written explicitly, are zero modes of the operators A, 
A + 3/C or (A + 3/C)^ and are therefore homogeneous modes or gauge modes which should 
be ignored. 

After these simplifications, we rewrite our partially-continued correlator as 

dp sinpx e-- 9pin^^^x') (41) 



c 87r^(p2 + 4) sinhx sinhpvr a_p 



and insert the expression (|30|) to obtain 

(42) 



dp sinpx e" 



■pn 



c 87r2(p^ + 4) sinhx sinhpvr 



gpiX)g^piX') + ^g^{X)g,{X') 

Oi—p 



We are now ready to perform the analytic continuation X = —i'^ — t under the integral. 
Under this substitution, we have 

g^,{X) ^ e^gi^{r) (43) 

where the Lorentzian Jost function (^p (t) is defined to be solution to the Lorentzian pertur- 
bation equation Ogpij) = {p^ + 4)(yfp (r) obeying (^^ (r) -^ e"*^"^ as r — >• — cx). Equation ( ^3[ ) 



follows by matching at large X. Like gp{X), the Lorentzian Jost function (^p (r) is analytic 
in the upper half p-plane. 
The correlator now reads 



I 



dp sinpx 



c 87r^(p^ + 4) sinhx sinhpvr 



e-^%^(r)A(r') + ^^p^(r)^p^(r') 

Qi—p 



(44) 



We want to distort the contour of integration C back to the real p axis. This is only possible 
provided g^p is analytic in the region through which the contour is moved. 

Analyticity of g-p in the strip < Q'(p) < 1 is proven as follows. One may re-express 
the Lorentzian perturbation equation Oq = {p^ + 4)g in terms of Lorentzian proper time 
t. The equation then takes the form t^q + A{t)tq + B{t)q = —p^q. The coefficients A[t) 
and B{t) have Taylor expansions in t about t = 0. One now seeks power series solutions 
q r-^ f^ {1 + qit + q2t^ + ...) and finds the appropriate indicial equation for s and recursion 
relation for the coefficients. Here we obtain s = ±ip, and a recursion relation that is non- 
singular as long as p is not an integer times i. In this situation, one is guaranteed that the 
series for q{t) converges within a circle extending to the nearest singularity of the differential 
equation [^, and even then the solution may be defined by analytic continuation around 
that singularity. In our case the first singularity occurs on the real axis when the background 
scalar field velocity (po is zero, when infiation is over and reheating has begun. Matching 
across that singularity is accomplished by switching from q to "^n, and will not present any 
difficulties. 



Since g-p is analytic in the desired region, we may now distort the countour C back to 
the real p axis, recovering the bound state term from the simple pole at p = iA. Once the 
integral is along the real axis we use symmetry under p — > —p to rewrite it as 



dp 



smpx 



-oo 167r^(p^ + 4) sinhx 



cothpvr \gi{r)g'lj,{r') + g'lp{r)g;;{r') 



9i{r)9Ur') 



9-pi 



r)9i{r'] 



+ 



1 



sinhpvr 



-9i{r)9i{r') 



a 



-p ^L 



a,- "■ - P^ > 



9%{T)gt 



■ (45) 



Note that for real p, g^^ij) is the complex conjugate of gpij') and a_p is the complex 
conjugate of a^, so the second term in square brackets is imaginary, but the first and third 
terms are real. 

We have the correlator in the Lorentzian region. Recall that the derivation assumed 
X — X' > 0. This translates to r' > r, and we have calculated the Feynman (time ordered) 
correlator for all x subject to this condition. For cosmological applications we are usually 
interested in the expectation value of some quantity squared, like the microwave background 
multipole moments or the Fourier modes of the density field. For this purpose, all that 
matters is the symmetrized correlator ({g(x, r), g(0, r')}) which is just the real part of the 
Feynman correlator. The symmetrised correlator also represents the 'classical' part of the 
two point function and in the situations of interest it will be much larger than the imaginary 
piece. 

The second term in (^5]) is pure imaginary and does not contribute to the symmetrized 
correlator. The other terms combine to give our final result 



({^^(X,r),vl/^(0,r')}) 
dp 



fi;20o(r)0o(r 
sinpx 



1 



1 



5R(cothpvr^^^(r)^^^(r') + ^.^j^^^ 



-^9'Ar)gl:{r' 



sinh Ax af^ ^f^ (r) ^f^ (^') 



47r(4-A2) sinhx M sinAvr 



(46) 



where the bound state pole arises as mentioned above from distorting the contour across the 
pole dX p = iK and we have converted the final integral along the real axis to one from to 
oo. We have also converted from q variable to \E'Ar by multiplying by the appropriate factor 
where dots denote derivatives with respect to Lorentzian proper time. Also recall that units 
here are such that the comoving curvature scale is unity (/C = —1). 

The first term in this formula is essentially identical to that derived by |^ and ||10|| , 
although in those derivations it was obtained only as an approximation. The second term, 
representing the reflection amplitude for waves incident from the regular pole, is only impor- 
tant at low p, since the sinhpvr denominator suppresses it exponentially at high p. Finally, 
the bound state term produces long range correlations beyond the curvature scale. For large 
p we recover the usual scale-invariant spectrum of inflationary quantum fluctuations. To see 
this note that at large p,\gp {t)\ ^ 1. Then according to (^61) there are equal contributions 
to the variance (\E'^(0, r)) from each logarithmic interval in p at large p. 

As mentioned in the introduction, one of our main concerns is to establish the differences 
between perturbations about singular and non-singular instantons. The above derivation was 
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for singular instantons. It is straightforward to follow the argument through for non-singular 
instantons, where the coordinate X now runs from — oo instead of zero. The only change 
in the final formula (|i6| ) is that the phase factor ap/a^p is replaced by Cp/dp, which is the 
reflection amplitude for waves incident on the potential from X = +00. 

VII. CONCLUSIONS 



We have derived a general formula p6|) for the time dependent correlator of the Newto- 
nian potential in open inflationary universes resulting from Euclidean cosmological instan- 
tons. The phase shifts and mode normalizations in the Euclidean region can be calculated 
numerically, and the Lorentzian Jost functions can be evolved numerically until they pass 
outside of the Hubble radius during inflation and freeze out. In future work we shall obtain 
the associated primordial matter power spectra and cosmic microwave anisotropies for a va- 
riety of potentials, considering both singular Hawking-Turok and non-singular Coleman-De 
Luccia instantons [^ . 



To summarise, we have computed from first principles the spectrum of density pertur- 
bations in an inflationary open universe including gravitational effects. We found that the 
Euclidean path integral coupled to the no boundary proposal gives a well defined, unique 
fluctuation spectrum, obtained via analytic continuation of the real-space Euclidean corre- 
lator. 
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